Electronic transport in quantum cascade structures 
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The transport in complex multiple quantum well heterostructures is theoretically described. The 
model is focused on quantum cascade detectors, which represent an exciting challenge due to the 
complexity of the structure containing 7 or 8 quantum wells of different widths. Electronic transport 
can be fully described without any adjustable parameter. Diffusion from one subband to another is 
calculated with a standard electron-optical phonon hamiltonian, and the electronic transport results 
from a parallel flow of electrons using all the possible paths through the different subbands. Finally, 
the resistance of such a complex device is given by a simple expression, with an excellent agreement 
with experimental results. This relation involves the sum of transitions rates between subbands, 
from one period of the device to the next one. This relation appears as an Einstein relation adapted 
to the case of complex multiple quantum structures. 



INTRODUCTION 



Quantum Well Infrared Photo-detectors (QWIPs) have 
become widely used quantum heterostructures for ther- 
mal imaging application during the last ten years 1 . 
Thanks to the very recent introduction of Quantum Cas- 
cade Detectors (QCDs)Sii4, QWIPs exist now in two 
modes of electronic transport : photoconductive and pho- 
tovoltaic. For the first type (the regular QWIP inserted 
today in focal plane arrays), an external bias is necessary 
to drift the electrons. They are characterized by a dark 
current and so they operate at low temperature. The 
second type of detectors uses a photovoltaic mode : no 
external bias is applied. No dark current exists in these 
devices and they are promising for applications at higher 
temperature, larger wavelengths 3 and smaller pixel areas. 
In all these structures, however, the electronic transport 
is a very complex issue and the performances of the de- 
vices are finally governed by fundamental physics prob- 
lems. 

In a standard QWIP, the transport involves both 2D 
and 3D electronic states in the quantum well and in the 
continuum, respectively. The diffusion of electrons from 
3D to 2D states (and vice versa) is a particularly dif- 
ficult theoretical problem^. That is why most mod- 
els use adjustable parameters such as the capture time? 
or the capture probability, affecting the photoconduc- 
tive gain 8 . In quantum cascade lasers, on another side, 
significant electric fields are applied, resulting in charge 
transfers. The description of the electronic transport 
has to involve Schrodinger and Poisson equations self- 
consistantly, which is a source of difficulty as well. With 
respect to these two devices, the QCD structure is rather 
simple since it involves only 2D states (and therefore ma- 
trix elements can be calculated without difficulty) and no 
bias is applied : there is no need for Poisson equation, 
and the flat band condition is an excellent approxima- 
tion, also because the doping levels are low. For all these 
reasons, the QCD appears as an archetype of complex 
quantum heterostructures, a model system for the ex- 
perimental study of the electronic transport in a com- 
plex multiple quantum well structure together with its 
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FIG. 1: Presentation of a period of a QCD. 1(a) Conduction 
band diagram and wave function associated with each energy 
level of a period. 1(b) Principle of a detection. 



modeling. This will be shown in this paper, in which the 
mechanisms of electronic transport in quantum cascade 
heterostructures are described. For the first time, elec- 
tronic transport in such a complex structure is calculated 
ab initio. By this expression, we mean that the current 
is calculated without any other adjustable parameters 
than standard effective masses and band offsets found 
in the general literature on semiconductors. The details 
of the QCD structure and the principle of detection are 
described in Sec. II. In Sec. Ill, the model describing the 
electronic transport is presented in details, relying on 
the calculation of all the transition rates between differ- 
ent subbands due to an electron-optical phonon hamilto- 
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nian. Through an analytical derivation of the sum of the 
different paths that an electron can follow to cross the 
structure, the resistance of the complex multiple quan- 
tum well structure is expressed in a very simple form 
(Sec. IV). The analogy between this expression and the 
Einstein relation will be underlined. 



II. PRESENTATION OF A QCD. 

The QCD considered in this paper is a GaAs/AlGaAs 
heterostructure composed of 40 periods of 7 quantum 
wells. The quantum wells are made in GaAs. The first 
quantum well of each period is n-doped in order to 
populate its first energy level E\ in the conduction band 
with electrons (the nominal doping concentration is 
about 5.10 11 cm -2 ). Barriers are made in Alo.34GaAso.66- 
The quantum wells (respectively the barriers) have the 
following widths: 68, 20, 23, 28, 34, 39 and 48 A 
(respectively 56.5, 39.55, 31, 31, 31, 31, 22.6 A). These 
dimensions have been optimized in order to create the 
quantum cascade of levels shown in figure 1. Figure 1 
(a) presents wave functions associated with each energy 
subband, in one period of the device. Wave functions 
associated with each energy subband are delocalized 
through the heterostructure. Absorption of a photon 
brings electrons from the first energy level to energy 
Ej or Eg- The first and the second quantum wells are 
coupled. In this way, electrons are transferred from level 
Ej or E$ towards energy level Eq mostly through the 
longitudinal optical phonon interaction. Because of the 
coupling between the second and the third quantum 
well, electrons will be relaxed to energy level E§ and so 
on. Electrons are relaxing from the left to the right of 
the structure step by step (see figure 1 (b)). Finally the 
net result under illumination is the transfer of electrons 
from energy level E\ of one period to energy level E\ of 
the following period. In this way, the structure acts as a 
photovoltaic infrared detector. 



III. MODEL OF THE ELECTRONIC 
TRANSPORT 

As usual with this kind of device, the resistance will be 
presented in terms of RqA (where Rq is the resistance of 
the pixel and A the area of the pixel) . Hypotheses about 
electronic transitions between the different energy levels 
have to be made for the determination of the RqA of such 
a structure. Considering the well and barrier widths, only 
interactions between electrons and optical phonons (LO- 
phonons) have been taken into account. The differences 
between the energy levels are indeed sufficiently high to 
neglect the influence of the interaction between electrons 
and acoustical phononsS. All other possible interactions 
have been neglected too. Interface roughness has low in- 
fluence on intersubband scattering, although it was pos- 



sible to measure its influence at very low temperature 10 
(4 K). Electron-electron interaction is efficient for elec- 
trons thermalization inside a subbandii*i£, but negligible 
for intersubband scattering in our cascade scheme with 
many subband separations greater than 36 meV which 
corresponds to the energy of the optical phononsi&ii. 
The transition rates due to the interaction between elec- 
trons and optical phonons will be evaluated first follow- 
ing Ferreira et alii^. Starting from an initial state of 
wavevector k and energy E in the subband i, the tran- 
sition rate S°'j e (E) towards the subband j (in s" 1 ) is ob- 
tained through the integration of a matrix element in- 
volving a standard electron-optical-phonon Hamiltonian. 
This integration involves all the possible final states of 
energy E ± Twjlo m the subband j, where lolo is the en- 
ergy of a LO-phonon, the plus or minus sign accounting 
for absorption or emission of LO-phonons, correspond- 
ing to superscript a or e, respectively. This evaluation 
is made in the parabolic approximation of the energy 
bands. Transitions rates Sf 4 and S% are linked by the 
following equation : 



S? j (E) = SUE + ?uJLo) 



(1) 



Finally, the global transition rate Gy between the 
subband i and subband j is the sum of the two tran- 
sition rates for absorption of LO-phonons (G?-), and 
emission of LO-phonons (G?-). In order to calculate the 
global transition rates G?- and Gfj, all the initial states 
of energy E are filled at thermal equilibrium by the 
Fermi-Dirac occupation factor /. Let us first consider 
the electronic promotion from a subband i to a higher 
subband j (i.e. j > i). The integration on all these 
states is now performed on the subband i: 



Gf, = 



S? 3 {E)f{E){l-f{E + huj LO )) 

ej -Hu> L o 

xn opt D(E)dE 



-Hex:- 



StAE)f(E)(l-f(E-nu LO )) 



x (1 + n opt ) D(E)dE 



(2) 



(3) 



^opt 



is the Bose-Einstein statistic function which ac- 



counts for phonon population, D(E) the two-dimensional 
density of states of the subband i and ej is the minimum 
of energy of the subband j. Of course, similar expressions 
can be written in the case of electronic transfers from a 
subband i to a lower subband j. Figure 2 summarizes 
the four possible electronic transitions. 



At thermodynamical equilibrium, each global coefficient 
Gij is equal to its reciprocal factor Gji. The number 
of electrons per second transferred from subband i to 
subband j is equal to the number of electrons per second 
transferred back from subband j to subband i. The net 
current is null and the population in every subband is 
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FIG. 2: The different types of electronic transitions between 
two subbands i and j 

stationnary. 

In the context of a photovoltaic detector, the resistance 
at V is a key parameter since it governs the Johnson 
noise of the device. The calculation of this resistance 
needs the determination of the current for a very small 
applied bias. The bias can then be introduced as a 
perturbation, and therefore we will suppose that the 
matrix elements and the S%j (E) transitions rates do 
not depend on the voltage. Indeed, the electronic 
wavefunctions given by the weakly perturbed potential 
can be supposed non dependant on the arbitrary small 
voltage. In the frame of this approximation, the current 
results only from the variation with the applied bias of 
the distribution of the population as a function of the 
energy which affect the GijS through the occupation 
factors /. 

As an illustration of the exchanges of electrons be- 
tween subbands occurring in the structure, Figure 
3 represents the main transition rates at thermo- 
dynamical equilibrium at 80 K. In addition, these 
transfer rates are indicated in table I. Transition times 
Tij can be deduced from the global transition rate Gij by: 

nj = ni/Gij (4) 

where m is the two dimensional carrier density associ- 
ated with the subband i. Two consecutive cascades A 
and B are represented in Figure 3. For sake of clarity, 
transitions between two levels in the same cascade have 
been represented separately from transitions between two 
levels in consecutive cascades. On Figure 3(a), most sig- 
nificant intra-cascade transitions are represented. The 
typical value of transitions rates between two neighbor- 
ing levels are between a few 10 20 and 10 25 m -2 s , corre- 
sponding to transition times between a few ps and a few 
tens of ps for a temperature of 80 K. On Figure 3(b), the 
inter-cascade transitions are reported. Again, the solid 
line arrows represent the main transfer rates between one 
cascade and the following one, but these transitions are 
now limited to a few 10 18 m^ 2 s _1 at the same tempera- 
ture (and a corresponding transition time greater than 1 
(is due also to the low amount of electron promotion to 
higher subbands for satisfying Fermi Dirac distribution). 



The dashed line arrows represent other minor transitions. 
The comparison between intra- and inter-cascade transi- 
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FIG. 3: Major transition rates of two consecutive cascades 
of the real device at 80 K. Figure 3(a) represents the main 
transition rates inside each cascade (only the transition rates 
greater than 10 20 m -2 S -1 are represented). Figure 3(b) rep- 
resents the main transitions between the cascades : solid lines 
concern the major transition (greater than 4. 10 18 m _2 s _1 ) 
whereas dashed lines represents the other main transitions 
(greater than 1. 10 18 m _2 s _1 and lower than the major tran- 
sitions) . 

tions shows that the electronic mobility is higher inside 
a cascade than between two consecutive cascades by sev- 
eral orders of magnitude. This justifies our classification 
between inter- and intra-cascade transitions and has ma- 
jor consequences on the understanding of the transport 
in the structure. It recalls the situation in a p-n junction, 
where intra-band thermalization is faster than interband 
recombination by several orders of magnitude (picosec- 
onds versus nanoseconds). More generally, this is a very 
general situation in the physics of transport where a sys- 
tem can be divided into two reservoirs separated by a bot- 
tleneck for the conduction. The concept of quasi Fermi 
level can be imported from the p-n junction formalism 
into our case. Each cascade stays at a quasi thermody- 
namical equilibrium. Two quasi Fermi levels associated 
with two consecutive cascades can be considered : the 
cascade A (resp. B) is associated with the quasi-Fermi 
energy Ep (resp. Ep). At thermodynamical equilibrium 
without applying any voltage these two quasi-Fermi lev- 
els are equal. The structure is characterized by only one 
Fermi level. In the presence of an external applied volt- 
age, a gap between the two quasi Fermi levels appears 
and increases with the voltage : Ep = Ep + qV. Each 
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TABLE I: Values of some transition rates in relation to the 
main electronic transition represented on the figure 3. The 
left column concern the higher transition rates inside each 
cascade (intra-cascade transition rate) and the right column 
gives the higher transition rates between two consecutive cas- 
cades (inter-cascade transitions rates). These values are eval- 
uated at 80 K. 



energy level of the cascade A increases of the same factor 
qV too. Figure 4 represents the situation of two consec- 
utive cascades under a bias voltage V. 
The main transitions at 80 K between two consecutive 
cascades occur between subband to the subbands 
Ef and E B (transition towards level Ef is more effi- 
cient than the other transitions). In conclusion, the re- 
sistance of a QCD is completely determined at 80 K by 
a few inter-cascade cross transitions, namely E^ — > E^ 
and Ei — > E B . The optimization of a QCD requires 
decreasing these transitions rates, thus increasing the re- 
sistance and decreasing the noise figure. This is possible 
by a separation of the wavefunctions, but at the expanse 
of a lower optical matrix clement, and a lower response. 
All the challenge of the QCD design consists of mastering 
this trade off. 



FIG. 4: Evolution of the first cascade in relation to the second 
cascade under a bias voltage. Solid lines represent the position 
of the energy levels under a bias voltage and dashed lines 
concern the position of the energy levels without any applied 
voltage. Fermi levels are represented on this figure in bold 
lines. 



Let us consider two subbands i and j associated respec- 
tively with cascades A and B. The introduction of two 
quasi Fermi levels implies the differentiation of the two 
Fermi occupation factors in Eq. [5] and Eq. 03 G% and 
can be evaluated by : 



Gfj(Y) = / S«(E)f A (E)(l-f B (E + hLo LO )) 

xn opt D{E)dE (6) 

r+oo 

G%(V) = / S] l (E)f B (E)(l-f A (E-hoj LO )) 

x(l + n opt )D(E)dE (7) 

where /a and f B are the Fermi-Dirac occupation factors 
associated with quasi Fermi level Ep and Ep . Consider- 
ing Eq. n the difference is then equal to: 



IV. DERIVATION OF AN EINSTEIN 
RELATION 



G%{V)-G%{V)- 



where: 



SUE)a(E)(l- 1 (E))dE 
(8) 



The global current density will be evaluated by 
counting the electronic transitions between the two 
consecutive cascades. On one hand, each cascade stays 
at thermodynamical equilibrium, but on the other hand 
equilibrium between the two consecutives cascade is 
broken. As a consequence, transition rates from cascade 
A to cascade B are not equal to the reciprocal transitions 
rates from cascade B to cascade A. A global current 
appears. This global current density is given by: 



(•5) 



a(E) = n opt f A {E){l- f B (E + hcu LO ))D{E) (9) 

(F) = f B (E + hu; L o)(l-f A (E))(l + n opt ) 
71 1 fA(E)(l-f B (E + nu> L o))n opt [ ' 

Without any applied voltage, the term a(E) is equal to 
a eq (E) given by: 



a eq (E) = n opt f(E) (1 - f(E + huj LO )) D(E) (11) 

It corresponds to the first term of the series expansion of 
a as a function of the voltage V. Expressing the Fermi 
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Dirac functions, j(E) is simplified into : 



where T is the temperature of the sample and the 
Boltzman constant. 

We recall that in the context of infrared photovoltaic 
detection, applied biases are very small. The Johnson 
noise is related to the resistance at OV: Rq. In this 
calculation, it is then justified to linearize : 

1-7(*Q«^V d3) 

This leads directly to the linear I(V) behavior of the 
structure at low bias, through the multiplication by the 
constant a eq (calculated with no applied voltage). For 
little variations of the voltage, the difference can be 
approximated by the following equation: 



Einstein relation. By an identification between RqA and 
I 2 / (gn2D/i) (where I is the period of the structure, n^D 
is the 2D-density and /j, is the electron mobility) we find 
indeed again D/fj, = k^T/q, with the diffusion coefficient 
D expressed classically as a function of the transfer rates 
at equilibrium Gij by : 




~i 1 1 1 r~ 

4 6 8 10 12 

1000/T(K ) 



Gtj(V)-G%(V) « / S^(E)a e "^=dE(U) 

qV 



£j-nu> LO 

G?JV = 0V) 



k h T 



(15) 



A similar expression can be found for the difference : 



G\ j {V)-G a ji {V)ttG%{V = 0V) 



qV_ 

hT 



(16) 



Finally, global current density is so evaluated by the 
formula: 



qV_ 

k b T 



(17) 



where the term Gij is defined by the sum of Gfj and 
G\j calculated without any applied voltage. RqA can be 
finally deduced from the last equation: 



RqA 



k h T 



-A />' ('Li 



(18) 



We have previously made the hypothesis that the 
electronic wavefunctions and the matrix elements are 
not dependant on the small bias. The current has been 
calculated as a result from a variation of the distribution 
of carriers as a function of the energy from cascade A 
to cascade B. This means that the current appears as a 
diffusion current, in a system with a non-homogenous 
chemical potential. In this context, it is expected to 
find an expression of the RqA which looks like an 



FIG. 5: RqA as a function of 1000/T, where T is the tempera- 
ture of the sample. Solid line corresponds to the experimental 
curve and the markers are the results of the presented model. 

We find again D = 1 2 /t used in the context of diffusion 
currents in p-n junctions^. It is remarkable in Eq. 1181 
that the RqA can be calculated in a structure with no 
applied bias, considering only the transition rates Gij at 
equilibrium, from one cascade to the next. Eq. 1181 ap- 
pears as the main result of this paper. We have shown 
that the electronic transport in a complex structure as a 
QCD can be simplified. The transport is entirely dom- 
inated by 2 or 3 matrix elements present in the sum 
in Eq. 1181 which is interpreted as an Einstein relation 
adapted to QCDs. 

Figure 5 presents a typical example chosen to illustrate 
the quality of the modeling : the experimental RqA of the 
device as a function of 1000/T, where T is the tempera- 
ture of the sample, and the result of our modeling using 
the previous Eq. Activation energy deduced from 

experimental curve lies around 120 meV. This energy cor- 
responds to the transition between the first energy level 
E\ and the energy level Eq. Dark current originates es- 
sentially from electrons transferred from level E\ to level 
Eq, which is in a good agreement with the evaluation 
of the transfer rates presented in this article. The theo- 
retical modeling provides a good fit of the experimental 
curve : activation energy deduced from the model lies 
around 120 meV too. 

The resistivity given by the model is slightly lower than 
the experimental value. We atribute this small discrep- 
ancy to a doping level slightly lower than expected. In- 
deed, including a lower doping level (around 3.10 11 cm -2 ) 
in the model provides a perfect fit of the experimental 
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result, and on the other hand the measurement of the in- 
tegrated absoprtion agrees very well with a lower value of 
the doping level with respect to the nominal value. Other 
small discrepancies between the experimental results and 
our model could come from other deviations of the struc- 
ture parameters from the nominal values. The influence 
of the variations of structure parameters on the resistivity 
will be systematically reported in a forthcoming study. 
This is indeed necessary to analyze the performances of 
QCD focal plane arrays, in term of homogeneity. 

V. CONCLUSION 

In conclusion, we have reported a model of electronic 
transfers in a multiple quantum well device. Lifetimes 



associated with all transitions between two subbands are 
evaluated by considering only interactions between elec- 
trons and optical phonons. The resistivity of such a struc- 
ture around a null bias voltage can then be deduced by 
Eq 1181 which appears as an Einstein relation. An ex- 
cellent agreement is found between experimental result 
and our model, which does not include any adjustable 
parameters. 
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